Mon. Not. R. Astron. Soc. 000.[flH2l(20071 Printed 3 February 2008 (MN BTgC style file v2.2) 



Inverse- Compton emission from the lobes of 3C 353 



J.L. Goodger, 1 * M.J. Hardcastle, 1 J.H. Croston, 1 N.E. Kassim 2 and R.A. Perley 3 

1 University of Hertfordshire, College Lane, Hatfield, Hertfordshire ALIO 9AB, UK 

2 US Naval Research Lab 4555 Overlook Ave., SW Washington, DC 20375, USA 
3 National Radio Astronomy Observatory, P.O. Box O, Socorro, NM 87801, USA 



ABSTRACT 

X-ray emission due to inverse-Compton scattering of microwave background photons by 
electrons in the lobes of powerful radio galaxies has now been seen in a large number of 
objects. Combining an inverse-Compton model for the lobe X-ray emission with information 
obtained from radio synchrotron emission provides a method of constraining the electron pop- 
ulation and magnetic field energy density, which cannot be accomplished using the radio data 
alone. Using six frequencies of new and archival radio data and new XMM-Newton observa- 
tions of the Fanaroff & Riley class II radio galaxy 3C 353, we show that inverse-Compton 
emission is detected in the radio lobes of this source at a level consistent with what is seen 
in other objects. We argue that variations in the X-ray/radio ratio in the brighter eastern lobe 
require positionally varying magnetic field strength. We also examine the X-ray nucleus and 
the cluster, Zw 1819.1-0108, spatially and spectrally. 
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1 INTRODUCTION 



Extended X-ray emis sion from the lobes of powerful, FRII 
jFanaroff & Rilevll 19741) radio galaxies and quasars is thought to 
be produced by inverse-Co mpton scattering of the cosmic mi- 
crowave background (CMB) i Feigelson et al. 19951: Tashiro et. all 



1 19981: iTashiro et alfeoOll; llsobe et al.ll2002klHardcastle et alj2002l) 
and o f infrared photons from the core dBrunetti. Setti & Comastril 
1 1997h . Although the synchrotron self-Compton emission from the 
lobes can also be modelled, it is negligible compared to the inverse- 
Compton e mission from other po pulations except in the very small- 
est lobes drlardcastle et alj |2002). Recent studies of the integrated 
X-ra y properties of radio lobes have been carried o ut for large sam- 
ples dCroston et"ai]|2005l ; iKataoka & Stawarj|2005l) . showing that 
the CMB is the dominant photon population in most cases. Com- 
bining an inverse-Compton/CMB model for the lobe X-ray emis- 
sion with information obtained from radio synchrotron emission 
provides a method of constraining the electron population and mag- 
netic field energy density, which cannot be accomplished using ra- 
dio data alone. Unlike the jets, the lobe material is not moving 
relativistically so there are no b eaming effects to consider during 
this analysis (e.g. Mackav 1973). Understanding the contributions 
of the electron densities of the synchrotron and inverse-Compton 
electron populations and the magnetic field allows us to estimate 
directly the total energy in the radio source, and thus the amount of 



E-mail: j.l.goodger@herts.ac.uk 



energy that can be transferred to the environment, without assum- 
ing equipartition, and also allows us to investigate the distribution 
of in ternal energy within the sourc e. 

lHardcastle & Crostonl (2005) recently carried out a spatially 
resolved X-ray inverse-Compton study of the bright X-ray and ra- 
dio source Pictor A. They used the variation of the X-ray/radio ra- 
tio across the lobes to investigate the variation in electron density 
and magnetic field strength throughout the source. They found that 
variations in either component alone could not explain the observed 
X-ray and radio properties of the lobe. Similar spatially resolved X- 
ray studies of other bright FR II radio galaxies are needed to follow 
up these results. Here we report on observations of the radio galaxy 
3C353. 

3C 353 is a FRII radio galaxy associated with the cluster 
Zw 1718.1-0108. Although it is one of the brightest extragalactic 
sources in the sky at low frequencies, it is relatively poorly studied, 
presumably d ue to its low declination and low Galactic latitude. 
ISwainl dl996t) observed 3C 353 at four radio frequencies with the 
NRAO Very Large Array (VLA) and published a study of the po- 
larisati on variation across the jets, w hich favoured a spine-sheath 
model (Swain, Bridle & Baumf fl998). Their VLA observations re- 
vealed filamentary structure within the lobes making 3C 353 an ex- 
cellent target for a spatially resolved study of electron distribution 
and magnetic fields. X-ray emiss ion from the c l uster associated 
with 3C353 was first detected by llwasawa et all j2000h in ASCA 
and ROSAT images. The radio source resides in a giant elliptical 
galaxy on the edge of the cluster, which Iwasawa et al. found to be 
bright in the X-ray with a luminous point source coinciding with 
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the radio galaxy's core. Iwasawa et al. determined global cluster 
temperatures and used optical observations to identify the cluster's 
members. They confirmed a redshift of z — 0.0304 for 3C 353, a 
redshift of z — 0.028 for three bright member galaxies, and used 
the velocity distribution to confirm 3C 353 as a member of the same 
system. 

In this paper we use new XMM-Newton and radio observations 
to investigate the nature of the electron distribution and the mag- 
netic fields in the lobes and hotspots of 3C 353. We also determine 
the cluster temperature, density and pressure profiles and exam- 
ine the cluster's interaction with the lobes of 3C 353. Throughout 
the paper we use a cosmology in which Ho = 70 km s _1 Mpc -1 , 
O m = 0.3, and Q\ — 0.7. The angular scale is 1 arcsec = 0.61 kpc. 
We define the spectral index in the sense that S v oc v~ a . The pho- 
ton index T = 1 + a. 



2 DATA 

We used a combination of new and archival radio data and new 
XMM-Newton X-ray data to examine the broad- spectrum electron 
energy distribution of 3C 353. Our radio data extends to lower fre- 
quencies than have been studied previously, whilst our X-ray data 
has greatly i mproved sensit i vity a nd resolution compared to the 
ASCA data of llwasawa etafl feOOOh . 

2.1 Radio 

We obtained VLA observations at 1.6 GHz and 4.8 GHz from the 
VLA archive as well as carrying out new observations at 327 MHz 
at the VLA and at 620 MHz at the NCRA Giant Meter- wave Radio 
Telescope (GMRT) in 2006. We also obtained images of 3C 353 
at 74 MHz and 8.4 GHz (see below). All the radio data were re- 
duced in AIPS. The radio observations are summarized in Table Q] 
and Fig.Q]shows an example of the image quality achieved. 

2.1.1 4 Band; 74 MHz, 4 m 

The 74-MH z obs ervations used are those discussed by 
Kass im et aD J2007b and the data reduction is described in 
that paper. 

2.1.2 P Band; 327 MHz, 91 cm 

Our 327 MHz observations of 3C 353 were taken at the VLA on 
the 27th May 2006 using the nearby Pie Town antenna of the Very 
Large Baseline Array (VLBA) to increase the long baseline of our 
A configuration data set by a factor of two. We carefully monitored 
the Pie Town antenna during the calibration steps to ensure that the 
long baselines survived. 

Initial bandpass calibration using 3C 286 resulted in smooth 
phase and amplitude variation across all channels. As a flux cali- 
brator, 3C286 is slightly resolved at this frequency. If the fractional 
bandwidth is great enough that the intrinsic v isibility of 3C28 6 
varies across the pass band, a model is required dLazio et alj2005h . 
At 327 MHz with the Pie Town antenna, we found that this effect 
is negligible so 3C286 was used as the flux calibrator without a 
model, effectively treated as a point source. 3C 286 was used to 
calibrate the phase across all baselines. The resultant calibrated 
data set was flagged to remove noise then spectrally averaged using 
SPLAT, averaging every 4 channels to reduce the number of chan- 
nels to 6, whilst producing a single source data set. The calibrated 



data included baselines up to 80 kA, twice what would have been 
achieved with only the A configuration. 

The B-configuration data set was bandpass calibrated using 
3C 286 to calibrate the channel gains, then flux using 3C 286 and 
phase calibrated with 1416 + 067. After flagging and phase cross- 
calibrating with the A configuration data set, the data sets were 
merged and a multi-facet deep cleaned map of 106 fields was pro- 
duced. The dynamic range of the final map is 580:1 with a resolu- 
tion of 6.7 x 2.9 arcsec. 



2.1.3 620MHz,48cm 

Our 620 MHz observations were taken at the GMRT on 19th July 
2006. The two observing frequencies (upper and lower sidebands) 
were calibrated separately. Preliminary bandpass calibration was 
performed using 3C 286 followed by flagging and one iteration of 
flgit to remove RFI from all sources. This flagged data set was then 
bandpass- and flux-calibrated before being manually flagged and 
another iteration of flgit applied. The data set was spectrally aver- 
aged into a multi-source data set, cutting the number of channels 
from 128 to 6, whilst applying the bandpass calibration. This data 
set was finally flux- and phase-calibrated before 3C 353 was split 
into a single-source file. The flux and phase calibrator used was 
1743-038. One iteration of phase self-calibration was performed 
before the data from the two observing frequencies were split into 
individual channels and recombined using DBCON to a single data 
set containing all of the averaged channels. This combined data set 
was phase self-calibrated twice more and the final image has reso- 
lution 7.50 x 4.41 and dynamic range 590 : 1. 

2.1.4 L Band; 1.67 GHz, 22 cm 

Observations from A-, B- and C-configurations were flux- and 
phase-calibrated separately using NRAO530 as the phase calibra- 
tor. Due to the allocation of the observing frequencies in the C- 
configuration data set, the data at the two frequencies were split 
and recombined using DBCON to match the frequency allocations 
of the A- and B-configuration data sets. Once combined, the data 
were phase self-calibrated and deep cleaned, with 100,000 itera- 
tions of CLEAN applied. The final image has a resolution of 1.7 x 
1.4 arcsec and a dynamic range of 500 : 1. 

2.1.5 CBand; 4.8GHz, 6.3cm 

For maximum uv coverage, we use B-, C- and D-configuration data 
sets, each of which were phase- and flux-calibrated and flagged. We 
combine two epochs of D-configuration observations to increase 
the signal to noise. 

The final combined data set was phase and amplitude self- 
calibrated, normalizing the gain. There was no significant reduc- 
tion in the maximum flux when compared to phase-only self- 
calibration, but the image quality was improved. The final image 
has aresolution 2.24 x 1.81 arcsec and adynamic range of 1180:1. 

2.1.6 X Band; 8. 4 GHz, 3. 6 cm 

The X band images of 3C 353 were kindly provided by Alan Bridle. 
The deep cleaned images have a dynamic range of 2500: 1 and a 
resolution of 1.30 x 1.30 arcsec. The details of the reduction are 
described bv lSwairJ jl996h . 
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Figure 1. 1.67 GHz total intensity image of 3C 353 at 1.8 arcsec resolution 



Table 1. Radio Observation Details 



Band 


Frequency 
(MHz) 


Telescope:Config. 


Date 


Duration 

(s) 


Bandwidth 
(MHz) 


Phase Calibrator 


Dynamic Range 


Res. 

(arcsec) 


4 


73.8 


VLA 


07/03/1998 


7740 


1.3 


3C405 


2090:1 


25.0 x 25.0 


P 


327.3 


VLA:A+PT 


27/05/2006 


25660 


3.125 


3C286 


580:1 


6.75 x 2.92 






VLA:B 


04/10/1998 


3010 


0.098 


1416+067 








614 


GMRT 


19/07/2006 


19847 


16.0 


1743-038 


590:1 


6.41 x 4.27 


L 


1665/1385 


VLAA 


19/05/1968 


20710 


12.5 


NRAO530 


500:1 


1.78 x 1.44 






VLA:B 


14/08/1968 


6640 


25.0 


NRAO530 










VLA:C 


16/09/1985 


1441 


50.0 


1741-038 






C 


4848/4898 


VLA:B 


14/08/1986 


6799 


50.0 


NRAO530 


1180:1 


2.24 x 1.81 






VLA:C 


22/08/1985 


1679 


50.0 


1741-038 










VLA:D 


15/11/1993 


532 


50.0 


1730-130 










VLA:D 


28/09/1984 


3360 


50.0 


1725+044 






X 


8440/8452 


VLA:BCnD 


13/03/1994 




12.5 




2500:1 


1.30 x 1.30 



2.2 X-ray 

We observed 3C 353 on the 25th August 2006 and 17th February 
2007 with XMM-Newton EPIC MOS 1 , MOS2 and pn cameras. The 
initial observation (Set 1) from 25th August 2006 yielded 39 522 s 
and 39 525 s for MOS 1 and MOS2 respectively and 34 042 s for the 
pn camera, whereas the second observation on the 17th February 
2007 (Set 2) was for 10473 s and 10487 s for MOS1 and MOS2 
and 5 653 s for pn. The pn camera was in Extended Full Frame 
mode for both observations and the MOS cameras were in Full 
Frame mode. 

The data sets were initially processed using SAS version 7.0.0. 
Standard EPIC MOS and PN pipelines and the standard filters 



#XMMEA_EA and #XMMEA_EP were applied. The MOS data sets 
were filtered to include single, double, triple and quadruple events 
(PATTERN^ 12) whereas the PN data set was filtered to include 
only single and double events (PATTERN ^4). Set 1 MOS data 
sets were free from flare events whilst the PN data set was filtered 
with a count threshold of 2 counts per second to remove a small 
flare event at the end of the observation. In Set 2, all cameras were 
affected by flaring; the MOS and PN data sets were filtered with 
count rate thresholds of 0.3 and 0.7 counts per second respectively. 
As only 1510 seconds remained of the PN data, this data set was 
excluded from further analysis. The net livetimes for Sets 1 and 2 
were 44 039 for MOS 1 and MOS2 and 26 052 for PN. 
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Figure 2. Radio contour images of 3C 353 at 327 MHz (top left), 620 MHz (top right), 1.67 GHz (bottom left) and 4.8 GHz (bottom right) all convolved to 
7.5x7.5 arcsec resolution. The contour levels are at 5crx(l,2,4,...)mjy/beam where 5a = 0.0106 for 327 MHz, 0.0113 for 620 MHz, 0.0045 for 1.67 GHz 
and 0.0017 for 4.8GHz 



The MOS and PN data sets were energy filtered to include 
energies between 0.5 and 5 keV before using the SAS command 
evselect to generate images from the filtered events files. The task 
eexpmap was then used to generate exposure maps without vi- 
gnetting. We chose not to include vignetting correction in these 
images as it incorrectly weights up the particle background at the 
edge of the field, creating artefacts in the image. The images from 
the separate detectors were normalized to the Set 1 PN count rate 
(in order to remove the PN chip gaps) and then combined. The ex- 
posure maps were combined then divided by the total filtered live- 
time. The combined normalised image was finally divided by the 
scaled combined exposure map. The combined image was Gaus- 
sian smoothed with a kernel of 10.4 arcsec to highlight the clusters 
and radio lobes using the CIAO command aconvolve. 

X-ray spectra were extracted using especget, a script that com- 
bines the evselect, arfgen and rmfgen commands in SAS, from re- 
gions defined to examine the emission from the northern and south- 
ern sub-cluster regions, and the east and west lobes. As the core and 
lobes are small regions near the pointing centre, so that vignetting 
is not important, we used local background subtraction. However, 
for the cluster, vignetting needs to be considered. The filtered event 
files were weighted using evigweight be fore spectra were ex tracted 
using a double subtraction method, (e.g. lArnaud et alj2002h . Back- 
ground te mplate files for the field of view were made using files 
created bv lRead&Ponmanl (2003) and scaled to the same particle 
background level as our observations. A local background region 
away from the galaxy and cluster was defined using DS9. A tem- 
plate background spectrum was subtracted from the source spec- 



trum events file to account for instrument and particle noise before 
a local background region was used to subtract residual background 
emission due to the differences in the Galactic/extragalactic back- 
ground levels of the source and scaled background datasets. Due 
to the low signal to noise of Set 2, only Set 1 could be used in the 
double subtraction method applied to the cluster. The spectra were 
then binned to 20 counts per channel after background subtraction, 
ignoring the first 20 channels for the MOS cameras and the first 50 
channels for the pn camera. 



3 RESULTS 

In this section we discuss the results of the X-ray spectral fitting 
and radio flux density measurements. X-ray fitting was carried out 
using XSPEC v 1 1.3 in the energy range 0.3 - 7.0 keV. Where ther- 
mal models were fitted we used a redshift of 0.03. 



3.1 Absorbing column density 

The c olumn density towards 3C 353 is uncertain. Ilwasawa et al] 
d2000h adopt a value of 1.0 x 10 21 cm based on the H I mea- 
surements of iDickev & Lockmanl dl990l) . However, they note that 
the visual extinction in the direction of 3C 353 would correspond 
to a column dens ity of a factor ~ 2 hig her. From the Galactic dust 
measurements of ISchlegel et al] i ll 9981) we estimate an Av of 1.4 
mag, which would correspond to a column density of ~ 2.6 x 10 21 
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Figure 3. Radio contours of 3C 353 at 1.67 GHz with a 0.3-7.0 keV XMM-Newton image (MOS 1, MOS 2 + pn) of the X-ray emission from the radio lobes 
and the cluster Zw 171 8. 1-0108 (left) Gaussian smoothed with a kernal of 10.4 arcsec; the same image zoomed in to show the inverse-Compton emission from 
the lobes (right), with the East radio/X-ray hotspot and West radio hotspot labelled. 



for standard Galactic gas/dust ratios. The true value is likely to lie 
somewhere between these extremes. 

We therefore used our new data to estimate a column den- 
sity directly. We extracted spectra for bright point sources around 
3C 353 using local background regions. In Set 1 there were 5 point 
sources bright enough for spectroscopy. We fitted power-law mod- 
els with free absorption to these sources and found a weighted mean 
N H = (1.64 ± 0.07) x 10 21 cm~ 2 . As this is consistent with the 
range estimated above, we adopt the mean value in the analysis that 
follows. 

The X-ray point sources are not shown in Fig.[3]and were ex- 
cluded from all spectra extracted where the cluster was examined. 
The emission from the radio galaxy was also excluded. 



3.2 The Nucleus 

In the X-ray, the nuclear region was defined with a radius of 40 arc- 
sec to include as many nuclear photons as possible while excluding 
the X-ray emission from the radio lobes. The background was es- 
timated from a local background region positioned above the radio 
galaxy so as not to include any emission related to the radio galaxy. 
In the radio, the nuclear region was much smaller, adjusted in each 
data set to be a close fit to the nuclear emission. The peak intensity 
was determined using jmfit for each frequency, the results of which 
are shown in Table|2] 

We initially attempted to fit single power-law and thermal 
models to this spectrum; however, the fits were poor (x 2 = 466.6 
for 259 d.o.f. and X 2 = 2242.5 for 259 d.o.f. r espectively). We 
theref ore fitted a model which has been shown jHardcastle et al.l 
2006) to provide a good fit to the nuclear X-ray emission from 



Table 2. Radio Core Fluxes 



Band 


Frequency 


Flux 


Luminosity 




(MHz) 


(Jy) 


(WHz" 1 sr" 1 ) 


P 


327.3/329.6 


0.097 


1.5 xlO 22 




620.5/633.3 


0.089 


1.4xl0 22 


L 


1665/1678 


0.116 


1.8xl0 22 


C 


4848/4898 


0.140 


2.2xl0 22 


X 


8440/8452 


0.149 


2.3 xlO 22 



narrow-line radio galaxies, consisting of a power-law component at 
Galactic absorption (see Section [3~7T l, a second power-law compo- 
nent with redshifted intrinsic absorption and a kT = 1 keV thermal 
component. The best fitting model had \ 2 = 268 for 237 d.o.f. The 
core spectrum is shown in Fig. [4] The photon indices of the unab- 
sorbed and absorbed power-law components were F = 1.49 ± 0.1 
and 1.33 ± 0.3 respectively, and the intrinsic column density was 
6.1 ± 0.6 x 10 22 cm -2 . The unabsorbed X-ray component has a 
1 keV luminosity density of 1.86 x 10 16 WHz -1 sr" 1 , whereas the 
absorbed X-ray component with the absorption removed and ex- 
trapolated to include energies from 2-10 keV has a luminosity of 
2.82 x 10 42 ergs -1 . We discuss the interpretation of the nuclear 
emission in Section |4~T1 

3.3 Lobe properties 

By considering the X-ray emission from regions encompassing the 
east and west lobes but excluding the hotspots and the core, we ex- 
tracted spectra from the combined X-ray data set using local back- 
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Figure 4. X-ray spectrum of the nucleus in the energy range 0.3-7.0 keV. 
The model plotted is the dual power-law + thermal model as described in 
Section l3~2l For clarity only the PN data are plotted. 



3.4 The nature of the hotspots 

Both the East and West hotspots are clearly detected in all radio 
frequencies. The West hotspot is undetected in the X-ray; how- 
ever, the East hotspot has an X-ray counterpart 6.5 ± 0.8 arcsec 
(4.0 ± 0.5 kpc) from the centre of the radio hotspot. This has too 
low signal to noise for spectral analysis to be performed, with only 
60 counts between 0.3 and 7.0 keV. In recent Chandra observations 
(Kataoka et al., private communication), the same feature is de- 
tected but the offset was measured as 4.6 ± 0.5 arcsec. A jet knot 
was also detected in the vicinity of the East X-ray hotspot suggest- 
ing that the difference in our offset measurements could be due to 
contamination from this knot. 

The inverse-Compton analysis performed on the lobes can- 
not be applied here as the electron populations responsible for any 
inverse-Compton X-ray and radio synchrotron emission reside in 
different regions. X-ray hotspots in radio galaxies of similar radio 
luminosity to 3C 353 are generally thought to be produced by syn- 
chrotron emission, and so this is likely to be the dominant emission 
process for this X-ray hotspot as well. The offsets measured be- 
tween the X-ray and ra dio hotspots are consistent with what is seen 
in som e other sources ( lHardcastle et al J [20021 |2007| ; lErlund et al.l 
l2007h . 



ground regions which we fitted with both thermal emission models 
and power laws. Assuming a circularly symmetric model for the 
core, we determined that ~ 6 per cent of the counts in the East 
lobe spectra are from the core whereas the core counts account for 
a quarter of those in the West lobe spectra. Using the XMM-Newton 
calibration files, we determined the fraction of the core spectrun ex- 
pected to be scattered into the West lobe region to be 4 per cent of 
the total core emission, and so we included a fixed component in 
the West lobe spectral fit consisting of the best-fitting core model 
with a normalisation fixed at 0.04 times that of the core. The details 
of the core spectrum are given in Section [3~!2l Due to the low signal 
to noise of Set 2, the West lobe was undetected in this data set and 
so our spectra were extracted from Set 1 only. 

The background emission from the cluster environment is the 
primary contributor to the uncertainty in the lobe flux densities. 
To best take account of this variation we used a truncated annular 
background region for the East lobe where the lobe appears to inter- 
act with the cluster. We investigated the effect of using local back- 
ground regions above and below the galaxy; the best fitting power- 
law models gave consistent photon indices but led to an increase 
of a factor 2 in the flux density of the East lobe. These background 
regions only take account of the emission from one of the two sub- 
clusters leading us to believe that a truncated annular background 
region is a more accurate measure of the cluster's contribution. For 
the West lobe, a local background region located to the south of the 
galaxy, at a similar distance from the brightest cluster emission took 
account of the emission from the southern sub-cluster sufficiently 
to allow a good fit to the spectra. 

The power-law model gave better fits than the thermal model 
for both lobes with x 2 = 39.6 for 34 d.o.f for the East lobe and 
X 2 = 23.1 for 24 d.o.f for the west lobe. The photon indices were 
T = 1.9 ± 0.4 and V = 1.2 ± 0.6 respectively, coresponding 
to a — 0.9 ± 0.4 and a = 0.2 ± 0.6. Figure [5] shows the best- 
fitting power-law model for the East and West lobe X-ray data. The 
low frequency X-ray spectral index, measured between 327 MHz 
and 1.67 GHz, is ~ 0.7 for both lobes, so that the X-ray spectral 
indices of both lobes are consistent with a CMB inverse-Compton 
scattering model. We discuss the nature of the lobes in Section l4~2l 



3.5 Cluster characterization 

Unfortunately, the southernmost section of the southern sub-cluster 
extends beyond the field of view of the X MM-Newton cameras ; 
however, a comparison to the ASCA image of llwasawa et al.l ( feOOOh 
shows that the peak of the X-ray surface brightness is included in 
our data. 

We extracted spectra from the northern and southern sub- 
clusters using a double subtraction method with annular local back- 
ground regions. The northern and southern sub-cluster regions are 
shown in Fig. [7] Thermal models fitted using mekal in XSPEC 
give acceptable fits (\ 2 = 258 for 175 d.o.f and x 2 = 544 
for 490 d.o.f for the northern and southern sub-clusters respec- 
tively) with kT — 3.3 ± 0.3 keV for the northern sub-cluster and 
kT — 4.0 ± 0.5 keV for the southern sub-cluster both fitted with an 
absorption of N h = 1-6 x 10 2 1 an" 2 . The best-fitting models are 
shown in Fig.l6l llwasawa et al.l J2000) determined the temperature 
of the combined sub-clusters to be kT — 4.3 ± 0.2 keV using a 
Galactic absorption of Nh = 1-0 x 10 21 cm" 2 whereas when the 
upper limit (1.6 x 10 21 cm~ 2 ) was used, the temperature dropped to 
3.94 ± 0.21 keV which is consistent with the mean temperature of 
the cluster. The difference in temperature of ~ 1 keV supports the 
idea that the two sub-clusters are originally separate components 
undergoing a merger. Taking a north-south slice through the cluster 
using rectangular regions of 208 x 104 arcsec from 17:20:50.074, 
-00:53:13.77 (Fig. [8j, there is no evidence for a shock and no in- 
crease in surface brightness, indicating a non-violent interaction be- 
tween these components. 

A thermal model was fitted to annular regions in the northern 
and southern sub-clusters. The results are shown in Table [3] The 
residual local background was accounted for using regions in a rel- 
atively source free area of the field of view. Using an annular back- 
ground region gave consistant temperatures for both sub-clusters. 
The northern sub-cluster shows a linear increase in temperature 
with radius with a best fit of T(r) = O.OlOr + 1.528 keV while 
the southern sub-cluster is isothermal within the errors. 3C353's 
East lobe lies at a radius of 250 - 350 arcsec from the centre of 
the northern sub-cluster but it may also be affected by the south- 
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Figure 5. X-ray spectra for the East (left) and West (right) lobes in the energy range 0.3-7.0 keV. The East lobe shows the MOS 1, MOS 2 and PN cameras 
with the power-law model while the West lobe shows only the PN data for clarity with the dual power-law + thermal contribution from the core as well as the 
best-fitting power-law model from the lobe emission as described in Section 1331 




Table 3. Temperature profiles for the northern and southern sub-clusters 





Northern 






Southern 




Radius 


kT 


X 2 /d.o.f. 


Radius 


kT 


X 2 /d.o.f. 


(arcsec) 


(keV) 




(arcsec) 


(keV) 




0-50 


1.8 ±0.3 


29/19 


0-75 


5.2 ±0.7 


537/348 


50-100 


2.2 ±0.3 


121/65 


75 - 150 


4.2 ±0.3 


638/444 


100-150 


2.8 ±0.4 


130/102 


150-200 


4.7 ±0.4 


250/167 


150-250 


3.5 ±0.3 


129/108 


200 - 250 


4.9 ±0.5 


296/183 


250 - 350 


3.5 ±0.5 


286/172 









ern sub-cluster. The annular temperature at this radius (kT — 
3.5 ± 0.5 keV) is consistant with the global northern sub-cluster 
temperature (kT = 3.3 ± 0.3 keV) as well as the temperature mea- 
sured by the slice at this radius (kT = 3.9±0.6 keV). We therefore 
adopt a temperature of 3.5±0.5 keV for the environment of 3C 353. 

As 3C 353 interacts mainly with the northern sub-cluster and 
as the southern sub-cluster is not fully covered within the XMM- 
Newton field of view, surface brightness and pressure profiles were 
constructed for the northern sub-cluster only. To determine the cen- 



tre of the northern sub-cluster we used a centroiding routine for 
each of the XMM-Newton cameras. Annular regions were used to 
extract a surface brightness profile, excluding the chip gaps, X-ray 
point sources, emission from 3C353 and from the southern sub- 
cluster and the position of the missing chip in the MOS 1 cam- 
era. An annular region was used to account for the residual local 
background emission. The surface brightness profiles for each cam- 
era were fitted with a model consisting of a f3 model convolved 
with the XMM-Newton PSF based on the on-axis parametrisation 
described in the XMM-Newton CCF files XRT1_XPSF_0006.CCF, 
XRT2_XPSF_0007.CCF and XRT3_XPSF_0007.CCF before being 
fitted with a j3 model (Fig.[9j. 

We examined the annuli in quadrants to highlight any den- 
sity variation in a north-south or east-west direction. The north- 
ern quadrants and the south-east quadrant were consistent. We ex- 
cluded the south-west quadrant from the final profile to ensure that 
X-ray emission from the radio galaxy, and any local disturbances in 
the cluster gas, were removed. The resulting surface brightness pro- 
files were consistent for all cameras. The joint best-fitting model, 
for all cameras, had a radius of 194 arcsec and /3 = 0.35. Consid- 
ering a lcr confidence interval for 2 interesting parameters, /3 was 
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Figure 7. Northern and southern sub-cluster regions from which spec- 
tra were extracted to determine the sub-cluster global temperatures. The 
crossed regions are those which were excluded from the analysis and in- 
clude the X-ray point sources and 3C 353. The crosses indicate the limits of 
the slice used. The inner edge of a large annular region masking the edge of 
the camera can be seen at the bottom of the image. 



unconstrained (though we consider values below 0.35 and above 
0.9 unrealistic), whereas the core radius of the northern sub-cluster 
was limited to between 184 and 412 arcsec. We discuss the results 
of this section in Section|4~3l 



4 DISCUSSION 

4.1 The nature of the core 

The double-peaked nuclear X-ray spectrum, consisting of a heav- 
ily absorbed component together with one with only Galactic 
absorption, is typical of what is observed in narrow-line radio 
galaxies (NLRG) (e.g. [Sambruna. Eracleous & M ushotzkv 1 19991; 
Grandi, Malaguti & Fiocchi 2006; Belsole . Worrall & Hardcastlel 
20061 ; iHardcastle. Evans & Crostonl2006b . However, when we plot 
the unabsorbed luminosity of the heavily absorbed nuclear com- 
ponent in 3C 353 against its 1 78-MHz total radio luminosity, as 
done bv lHardcastle et all J2006b (their figure 3) we find that 3C 353 
lies almost 1 order of magnitude below the correlation they de- 
termined for NLRGs: that is, 3C 353 is overluminous in the radio 
for its X-ray luminosity or undeiiuminous in the X-ray for its ra- 
dio luminosity, and lies closer to the region of the plot occupied 
by low-excitation radio galaxies (LERGs). By contrast, the unab- 
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Figure 8. Temperature variation of north-south slice through cluster 
with rectangular regions of 208 X 104 arcsec from a zero of offset at 
17:20:50.074, -00:53:13.77. 
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Figure 9. Surface brightness profile for the northern sub-cluster showing 
the PN data set with the best- fitting fi model which has j3 = 0.35 and core 
radius 194 arcsec. 



sorbed component is consistent with the correlation found by Hard- 
castle et al. between unabsorbed X-ray luminosity and nuclear 5- 
GHz radio emission (their figure 3). 

The optical emissi on-line classification of 3C353 is uncer- 
tain. lLaing et al] d 1994b use the ratio of the flux in the [O III] 
line to that in the Ha to distinguish between the NLRGs and 
LERG s. The best pub lished optical spectrum of 3C 353 is that of 
Simps on" et alj d 1996b . and their quoted line fluxes clearly place 
3C353 below the [Oiii]/Hq cutoff of 0.2 proposed by Laing et 
al: thus by this definition 3C 353 would be a LERG. However, if 
the line fluxes are corrected for reddening us ing values determined 
from the Schlegel. Finkbei ner & Davisl i 19981) dust maps (see Sec- 
tion 13. IK then the emission-line ratio rises to 0.17 ± 0.02, where 
the errors are derived from the normal 10 per cent errors of Simp- 
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son et al. and are almost certainly too low: within the errors we 
cannot say whether 3C 353 should be classed as a LERG or NLRG. 
The nature of 3C 353's nucleus both in the optical and in the X-ray 
remains ambiguous. 



4.2 Lobe emission and properties 

Using the method of ICroston et al.l j2005t) . the X-ray observations 
were compared to the predictions of an inverse-Compton model, 
based on synchrotron m odelling of the radio data sets using SYNCH 
(Hardcastle et al. 1998). The broad-band spectra for the lobes are 
shown in Fig. [10] The measured 1-keV flux densities of the East 
and West lobes are 11.6 ± 1.6 and 2.7 ± 1.0 njy. The equipartition 
magnetic field strengths for the East and West lobes are 0.89 nT 
and 0.84 nT respectively, while if we assume that all the X-ray 
emission is inverse-Compton in origin then the measured flux den- 
sities correspond to magnetic field strengths of 0.44 ± 0.06 nT 
and 0.68 ± 0.2 nT. Therefore, while the weakly detected West 
lobe has a measured magnetic field strength within errors of the 
predicted equipartition value, the East lobe's measured magnetic 
field strength exceeds the predicted equipartition value by a factor 
~ 2 so that B ts/B eq ~ 0.4. The West lobe's measured mag- 
netic field strength is also consistent with a substantial departure 
from equipartition. This implies that the lobes of 3C 353 are elec- 
tron dominated, which is not unusual considerin g the range of field 
strengths in the sample o f ICroston et all [2005). This factor would 
be reduced for 3C 353 if the lobes were not in the plane of the sky. 

The measured magnetic field strengths in the lobe correspond 
to internal pressures of 3.4 ± 0.6 x 10~ 13 Pa and 1.9 ± 0.2 x 
10 -13 Pa for the East and West lobes respectively. (If the alternative 
local background regions discussed in Section [331 were used, the 
magnetic field strengths would be reduced by ~ 0.1 nT and the 
internal pressure of the East lobe would double, but this does not 
affect our conclusions.) 

Using the radio data sets, we constructed a spectral index map 
between the 1.67-GHz and 327 -MHz images, FigQT] It revealed a 
variation of Aa ~ 0.2 across both lobes. Within the radio lumi- 
nous region of the East lobe, excluding the hotspots, the spectral 
index is roughly constant, a = 0.66 ± 0.01 despite the filamentary 
structure seen in Fig.[TJ whereas in the West lobe, a = 0.63 ± 0.01 
excluding the hotspot. We note that the radio lobes do not appear 
to be entirely separate, and thus consider the region in between the 
radio luminous lobes, north and south of the core to be an inter- 
lobe region, which cannot be unambiguously associated with ei- 
ther lobe. This inter-lobe region has a relatively steep spectrum of 
a — 0.75 ±0.01 whereas the hotspot in both lobes exhibits a flatter 
spectrum of a — 0.60±0.01. Fig.ll llshows the spectral index map 
between 1.67 GHz and 327 MHz with contours from the 327-MHz 
map. 

The X-ray/radio ratio was determined for the hotspot, lobe 
and inter-lobe regions, using the high signal to noise of the X-ray 
data set. The 327 MHz radio map of resolution 7.0 x 7.0 arcsec 
was used with the convolved, exposure corrected X-ray map of 
energies 0.3 to 7.0 keV. The X-ray/radio ratio was found to be a 
factor of 4 greater in the steep inter-lobe region than in the flat, 
hotspot regions. If the magnetic field strength and the number den- 
sities of both the inverse-Compton and synchrotron emitting elec- 
tron are constant across the lobe, the X-ray/radio ratio would also 
be constant, but this is not what we see in 3C 353. Si milar results 
have been seen in oth er radio galaxies, notably 3C 452 dlsobe et all 
|2002|) and Pictor A (Har dcastle & Crostonl l2005h . Hardcastle & 



Croston consider the following three models to explain the vary- 
ing X-ray/radio ratio: 

(i) Some other emission process could boost the X-ray emission 
in the inner regions; 

(ii) The central regions contain more low-energy electrons rela- 
tive to the outer regions; 

(iii) The magnetic field strength varies as a function of position. 

A contribution from the core was included in the X-ray spec- 
tral analysis which also takes account of the contribution of the un- 
absorbed X-ray emission associated with the jet and the absorbed 
X-ray emission from the accretion disk. Emission from the core 
cannot therefore boost the X-ray emission in the inner regions. 
Our spectral fits show no evidence for a contamination of the lobe 
spectra by galaxy-scale thermal emission. Neither could a boost be 
due to the galaxy's proximity to the cluster as local background 
regions were used to account for any thermal emission from the 
sub-clusters. For inverse-Compton scattering of nuclear photons to 
be the dominant process, the core bolometric luminosity needs to 
be at least 10 38 W. This is not unrealistic, however mod elling the 
lobe su rface brightness for this scenario using th e results o flBrunettil 
d2000h as described bv lHardcastle et all d2002l) and with 3C 353 in 
the plane of the sky, reveals a gradient across the lobe which is not 
observed in either the X-ray image nor the X-ray /radi o ratio. Using 
the spectral energy distributions of lHaas et al.l d2004l) . we deduced 
a typical IR spectrum for 3C 353 by scaling a dual power-law fit 
to the spectral energy distribution of 3C 33. 3C 33 was chosen as it 
has a similar luminosity and is classified as a NLRG/FR II. The IR 
spectrum was normalised by the ratio of low-frequency radio lumi- 
nosities of the two sources. The predicted flux density for the nu- 
clear inverse-Compton emission is 6.5 x 10 -12 Jy at 1 keV, which is 
a factor of ~ 700 fainter than the predicted CMB inverse-Compton 
emission and the observed flux density (4.4 and 11.6 njy respec- 
tively). Thus we can rule out model (i). 

If we consider the magnetic field to be constant at the mea- 
sured lobe averaged magnetic field strength of 0.39 nT and apply 
an inverse-Compton model to the lobes, we find the emission at 
327 MHz traces electrons with 7 ~ 4000 whilst the measured 
inverse-Compton emission traces electrons at 7 ~ 1000. As the 
critical frequency for synchrotron emission goes as 7 2 , we find that 
a variation in the X-ray/radio ratio of a factor of 4 requires a varia- 
tion in the spectral index between 10 MHz and 327 MHz of ~ 0.5. 
Even if the equipartition magnetic field strength of 0.89 nT is as- 
sumed so that the 327 MHz emission traces electrons of 7 ~ 3000, 
the observed spectral index variation requires a factor ~ 1.2 vari- 
ation in the X-ray /radio ratio which is still much lower than ob- 
serve d. This is also the case with Pictor A dHardcastle & Crostonl 
2005); a variation in the low-energy electron densities alone (model 
ii) cannot explain the variation of X-ray/radio ratio across the lobe. 

Alternatively, if we assume constant electron densities for both 
the synchrotron and inverse-Compton emission electrons and also 
that the magnetic field strength does not vary along the line of 
sight, we find that the observed variation in the X-ray/radio ratio 
then requires a variation in the magnetic field of at most a factor of 
2.5. For a given frequency, this means the spectral index observed 
at low frequencies, between 1.67 GHz and 327 MHz, for regions 
of high X-ray/radio ratio should correspond to the spectral index 
between 3.5 GHz and 825 MHz for low X-ray/radio ratio regions. 
We consider the spectral indices between 4.8 GHz and 1.67 GHz to 
limit the spectral indices of the low X-ray/radio ratio regions and 
find that they exceed the upper limit predicted by the observed X- 
ray/radio ratio. From this comparison, we cannot rule out the pos- 
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Figure 10. Broad-band spectrum for the East (left) and West (right) lobes of 3C 353 . A synchrotron emission model (solid line) is fitted to the radio data 
(asterisks). The measured X-ray flux is represented by the dot with the bowtie indicating the error. The dotted line shows the inverse-Compton X-ray flux 
prediction and the dash-dotted line shows the predicted synchrotron self-Compton emission, both with the equipartition magnetic field strength. 



sibility that a varying magnetic field alone could be responsible for 
the observed variations in the X-ray/radio ratio and the spectral in- 
dices. We therefore considered a colour-colour d iagram using the 
method o f lKatz -Stone, Rudnick & Andersonl ( ll993h . In their study 
of Cygnus A, they argue that a single curve on a colour-colour 
diagram is consistent with a homogeneous distribution of rela- 
tivistic electrons together with varying the magnetic field strength 
across the source. The colour-colour diagram for 3C 353, shown in 
Fig. [12] has a simila r single curve and so is con sistent with this 
picture. In Pictor A, Hardcastle & Croston (2005) argued that the 
detailed positional variations of radio spectra and radio/X-ray ratio 
were not consistent with such a model, but our data are not good 
enough to rule it out here. 

We conclude that a varying electron spectrum alone cannot 
account for the observed variation in X-ray/radio ratio in the lobes 
of 3C 353 but that a magnetic field strength that varies by a factor of 
~ 2.5 throughout the lobes can explain it. Similar conclusions were 
reached by Hardcastle & Croston ( 2005) in their study of Pictor A. 
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Figure 11. Spectral index map between 1.67 GHz and 327 MHz with radio 
contours at 327MHz for levels of 0.01x(l,2,4,...)mjy/beam. The beam 
size is shown by the circle in the bottom left corner and the colour bar at the 
top shows the mapping of the grey levels to spectral index. 



4.3 The cluster-lobe interaction 

The surface brightn ess profile was converted to a pressure profile 
using the method of Birkins haw & Worralll dl993l) . so that the radio 
lobe pressures determined with SYNCH could be directly compared 
to the external pressure from the cluster at the position of the ra- 
dio galaxy. The centre of the East lobe lies ~ 280 arcsec from the 
centre of the northern sub-cluster. At this radial distance where the 
temperature of the environment is taken as 3.5 ± 0.5 keV, the exter- 
nal pressure is much greater than the internal pressure of the lobes 
(See Fig. 1131). As FRII radio galaxies are expected to be in pres- 
sure balance or over pressured, either 3C 353 is not in the plane 
of the cluster or there is an additional contribution to the pres- 
sure from non-radiating particles such as hot thermal or relativis- 
tic protons in the lobes. Previous studies of FRII radio galaxies in 
groups and clusters have shown that pressure balance can usually 
be achieved without additional protons (e .g. ICroston et al.l 120051 : 
lHardcastle e7al]|2002l ; iBelsole et al]|20"04r) . We therefore assume 



that 3C 353 is not in the plane of the cluster and using the East 
lobe, determine that 3C 353 requires a shift in radial position corre- 
sponding to ~ 600 arcsecs to be in pressure balance with the north- 
ern sub-cluster (here we use the best estimate of the lobe pressure 
from Section l4!2t . This places 3C 353 ~ 370 kpc either in front or 
behind the centre of the cluster. At this position, the West lobe's in- 
ternal pressure is also consistent with the intracluster pressure. The 
position of 3C 353 with respect to the pressure profile is shown 
in Fig.[T3] At this distance from the cluster, the external pressures 
seen by the lobes are different by a factor of ~ 1.8, which may help 
to explain the difference in appearance of the lobes, the East hav- 
ing a spherical appearance whilst the West lobe is elongated with 
hotspots at the outer edge. 
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Figure 12. Colour-colour diagram for 3C 353 made at 7 arcsec resolution. 
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Figure 13. Pressure profile for the northern sub-cluster shown by the greyed 
in region, with the measured pressure of the lobes of 3C 353 at the projected 
distance (solid line) and at the position with the radial shift applied to where 
the East lobe is in pressure balance with the sub-cluster (dashed line). The 
measured inverse-Compton emission was used to determine the pressure 
within each lobe. The temperature of the cluster at the radial distance of 
3C 353 was used. 



pendent of frequency we can rule out an foreground absorber and 
the required geometry is unrealistic for an obstruction in the lobe. 
With out further informatio n we are unable to identify this feature. 

Ilwasawa et alj ( feOOOh included optical observations of the 
cluster region centred on 3C 353 taken with the University of 
Hawaii 2.2m telescope in their analysis of Zw 1718.1-0108. They 
identified three additional massive galaxies, none of which reside in 
the northern sub-cluster. Deeper optical observations would help to 
establish whether there are any galaxies associated with the north- 
ern sub-cluster and whether 3C 353 is the dominant member of this 
sub-cluster despite its position at the edge. 



5 SUMMARY AND CONCLUSIONS 

Our results can be summarized as follows: 

• By fitting an inverse-Compton model to the lobes, we found 
the East lobe to be electron dominated and the West lobe to be 
consistent (within the large errors) with equipartition. 

• We determined that a variation in the electron spectron cannot 
account for the varying X-ray and radio emission alone, but that a 
change in the magnetic field strength across the lobes is required. 

• We have obtained a good X-ray spectrum of the nucleus of 
3C 353. Both the X-ray and optical properties of this source are 
ambiguous but it appears to lie in a region of parameter space in be- 
tween those normally occupied by narrow-line and low-excitation 
radio galaxies. 

• We have detected an X-ray counterpart for the East hotspot 
offset by 4.0 ± 0.5 kpc. 

• The northern and southern sub-clusters were found to be 
isothermal with a temperature difference of ~ IkeV supporting 
a model in which they are two originally separate components un- 
dergoing a merger with no evidence for a violent interaction. 
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4.4 Additional features of note 

In addition to the usual radio-lobe structure, the East lobe of 3C 353 
contains a dark circular region of unknown origin at 17:20:36.088, 
-00:58:44.15 (Fig. [TJ. Detectable in all radio frequencies except 
74 MHz (presumably due to its low resolution), the measured fluxes 
in this region are a factor of 100 above the rms background but are 
a factor of 2 fainter than the surrounding lobe emission. The X-ray 
image shows no sign of any feature in this region and neither does 
the Digital Sky Survey (DSS). As the level of the deficit is inde- 
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